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Abstract 



A simple confining separable interaction Ansatz for the rainbow-ladder 
truncated QCD Dyson- Schwinger equations is used to study quark deconfine- 
ment and meson states at finite temperature. The model is fixed at T = to 
implement quark confinement while preserving the Goldstone mechanism for 
the 7r. Within the Matsubara formalism, a very slow temperature dependence 
is found for the ir and p meson masses Mjj(T) until near the deconfinement 
£-~. . temperature T c = 143 MeV. Related to rapid decrease of the dynamically- 

generated quark mass function for T > T c , this model produces ir and p 
masses that rise significantly and are better interpreted as spatial screening 
masses. The T-dependent screening mass defect AMh = 2irT — Mh{T) is 
compared to results of lattice gauge theory simulations and also to those of 

^ an infrared dominant analytic model. 
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£2 1 Introduction 



The experimental search for the QCD deconfinement phase transition in ultrarel- 
ativistic heavy-ion collisions will enter a new stage when the relativistic heavy-ion 
collider (RHIC) at Brookhaven will provide data complementary to results from the 
CERN SPS Q. It is desirable to have a continuum field-theoretical modeling of 
quark deconfinement and chiral restoration at finite temperature and density (or 
chemical potential [/,) that can be extended also to hadronic observables in a rapid 
and transparent way. Significant steps in this direction have recently been taken 
through a continuum approach to QCDj^ based on the truncated Dyson-Schwinger 
equations (DSEs) within the Matsubara formalism M. For a recent review 

see ||. A most appealing feature of this approach to modeling nonperturbative 
QCDj- iA1 is that dynamical chiral symmetry breaking and confinement is embod- 
ied in the the model gluon 2-point function constrained by chiral observables at 
T = fi = and no new parameters are needed for extension to T, /J, > 0. Approx- 
imations introduced by a specific truncation scheme for the set of DSEs can be 
systematically relaxed. However due partly to the discrete Matsubara modes, the 
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finite T, /i extension of current realistic DSE models entails a complicated set of 
coupled integral equations. In the separable model we study here, detailed realism 
is sacrificed in the hope that a few dominant and essential features may be captured 
in a simple and transparent format. We simplifiy an existing T — \i — confining 
separable interaction Ansatz || by using a gaussian form factor. A related gaussian 
separable model has recently been explored for meson properties 0. 



2 Model Dyson-Schwinger equations at T = 

In a Feynman-like gauge where we take D^ v = 5 flu D(p — q) to be an effective in- 
teraction between quark colored vector currents, the rainbow approximation to the 
DSE for the quark propagator S(p) = [i$A(p) + B{p) + m ] _1 yields in Euclidean 
metric 



r> ( \ 16 f A n/ v B{q) + m 

B(P) = YJj2^ D{p - q) 



(2tt) 4 ^ ^qiA 2 (q) + [B(q)+m } 
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We study a separable interaction given by 

D(p -q)=D h{p 2 )h{q 2 ) + Di fi(p 2 )(P ■ q)h{<?) , (3) 

where D , Di are strength parameters and the form factors, for simplicity, are here 
taken to be fi(p 2 ) = exp(— p 2 /A 2 ) with range parameters Aj. It is easily veri- 
fied that if Do is non-zero, then B(p) = Am fo(p 2 ), and if D\ is non-zero, then 
A(p) = 1 + Aa fi{p 2 )- The DSE then reduces to nonlinear equations for the con- 
stants Am and Aa. The form factors should be chosen to simulate the p 2 dependence 
of A(p) and B(p) from a more realistic interaction. We restrict our considerations 
here to the rank-1 case where D\ = and A(p) = 1. The parameters Do, Ao and mo 
are used to produce reasonable n and uj properties as well as to ensure the produced 
B(p) has a reasonable strength with a range A ~ 0.6 ... 0.8 GeV to be realistic 0. 

If there are no solutions to p 2 A 2 (p) + (B(p) + m ) 2 = for real p 2 then the 
quarks are confined. If in the chiral limit (mo = 0) there is a nontrivial solution 
for B(p), then chiral symmetry is dynamical broken. Both phenomena are present 
in the separable model. In the chiral limit, the model is confining if Do is strong 
enough to make Am/Ao > l/y/2e. Thus for a typical range Ao, confinement will 
typically occur with M(p w 0) > 300 MeV. 

Mesons as qq bound states are described by the Bethe-Salpeter equation which 
in the ladder approximation for the present approach is 

-X H (P 2 )T H (p,P) = y^[D(p-q) 1 ,S(q + j)T H (q,P)S(q-^) 1 ,] , (4) 
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where P is the total momentum and H denotes a particular hadron. The quantity 
Xh is an eigenvalue and at the physical mass M# it satisfies Xh(P 2 = —Mh) = 1- 
With the rank-1 separable interaction, only the 75 and the 75P covariants contribute 
to the 7r H, and here we retain only the dominant form T n (p, P) = i / y 5 E n (p, P). For 
the vector meson, the only surviving form is T pp (p, P) = •y^(P)E p (p, P), with tJ(-P) 
being the projection of 7 M transverse to P. The separable solutions have the form 

. The 



Ei(p,P) = fo( P 2 )Q(P 2 



Ti, p , where the Ci factor out from Eq. 



eigenvalues A#(P 2 ) are the only determined quantities, and for example, 



K{P 2 



lQDr 



d A q 

W) 1 



p2 

(<? 2 - j)4 ff v + a s a s 



(5) 



The quark propagator amplitudes are defined by S(p) = —ip'crv(p 2 )+< J s(p 2 ), an d 
the d= superscripts indicate momentum arguments q ± P/2 respectively. For the 
transverse and longitudinal p states, there are expressions for Aj(P 2 ) and A^(P 2 ) 
analogous to Eq. (|5|). 

With parameters m /A = 0.0096, D Al = 128 and A = 0.687 GeV, this model 
yields M n = 0.14 GeV, M p = M u = 0.783 GeV, f n = 0.104 GeV, a quark condensate 
(qq) 1 / 3 = —0.248 GeV, and a p — 7 coupling constant g p = 5.04. 

In the limit where a zero momentum range for the interaction is simulated by 
foil 2 ) ^ 4 (<?)) then the expressions for the various BSE eigenvalues Xh(P 2 ) re- 
duce to those of the Munczek and Nemirovsky |§ model which implements extreme 
infrared dominance via 



D{p-q) = y-(27i) 4 ri 2 5^\p-q) . 
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(6) 



The correspondence is not complete because the quark DSE solution in this model 
has A{p) 7^ 1. We use r\ = 1.107 GeV to produce the same M u . The finite tempera- 
ture and chemical potential generalization of this infrared dominant (ID) model has 
been studied in [[|. 0- 



3 Finite temperature extension 

The generalisation to T 7^ is systematically accomplished by transcribing the 
Euclidean quark momentum via q — > (uj n ,q), where u n = (2n + 1)txT are the dis- 
crete Matsubara frequencies. Thus the integration measure in previous expressions 
becomes 




At finite temperature, the 0(4) symmetry is broken and the employed mass shell 
condition must be specified. One expects that if there is a bound qq meson state, 
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then the associated pole contribution to the relevant propagator or correlator will 
have a denominator proportional to 

1 - \ H (Q 2 m , P 2 ) ex Q 2 m + P 2 + M 2 H . (8) 

We investigate the meson mode eigenvalues Xh using only the lowest meson Matsub- 
ara mode (fl m = 0) and the continuation P 2 — > —M\. The masses so identified are 
spatial screening masses corresponding to a behavior exp(— Mux) in the conjugate 
3-space coordinate x and should correspond to the lowest bound state if one exists. 
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Figure 1: Spatial meson masses of tt and p for the Gaussian separable model (left 
panel) and for the analytic ID model (right panel). 

The results are displayed in the left panel of Fig. 1 along with the strength 
(Am) of the dynamical mass function at p — 0. The deconfinement temperature is 
indicated. The ir and p masses are only weakly dependent upon T until near the 
transition. For the tt this agrees with similar behavior obtained in the DSE approach 
with a more realistic interaction ||. For both the 3-space transverse and 3-space 
longitudinal p, the separate but still weak temperature dependence agrees with the 
findings 0] from the closely related ID model. The right panel of Fig. 1 displays the 
7i and transverse p results from the ID model in the chiral limit. In both panels of 
Fig. 1 we display the T-dependence obtained from the eigenmass condition Xh = I 
for T > T c in comparison with 2nT. 

We interpret the model results rising similar to 2nT in the following way. In the 
chiral limit and just above the transition, the dynamically generated mass function 
of the quarks has disappeared and the axial vector Ward identity that ties the pion 
Bethe-Salpeter amplitude in the 75 channel to this quantity indicates that there 
should be no solution of the pion BSE there. However the particular model used 
here to estimate Xh assumes the existence of a bound state with an amplitude 
or wave function profile given by the chosen form factor fo(q 2 ). The form of the 



4 



latter was chosen at T = to simulate the qualitative features of the DSE (and 
hence 7r BSE) solution. If bound state conditions disappear with rising temperature 
as interactions weaken, this approach is not well suited to detect it. The high 
temperature form of the denominators of the quark propagators in the loop integral 
for the "polarization" function A# is io 2 n + (q ± P/2) 2 . The dependence of the form 
factor on u 2 emphasizes the lowest mode. For ojq » A the form factor dictates 
that non-zero relative momentum q makes a minor contribution to the location 
of the lowest singularity in P 2 that determines the spatial screening mass. The 
higher the temperature the more the model behaves as if the fixed 3-space profile 
of the amplitudes had support effectively concentrated at q = 0, that is, almost 
independent of relative separation of quark and antiquark. Thus an approach to a 
screening mass characteristic of an independent fermion pair is to be expected here. 
In Fig. 2 this point is emphasized via comparison with spatial screening masses taken 
from || and representing lattice simulations above the transition. The horizontal 
lines mark 2tt, which in the right panel is corrected || for the lattice time extent 
N t . 
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Figure 2: Hadronic screening masses from the gaussian separable and ID models 
(left panel) compared to those of lattice gauge theory simulations (right panel, data 
taken from |§) 



From Figs. 1 and 2 it is evident that for T > T c , the "screening mass defect" 
AMh =2ttT — Mh{T) is small and decreasing rapidly in the separable model, it is 
apparently larger in the lattice simulations especially for the n, and it is significantly 
stronger with a larger temperature range in the ID model. The latter we attribute 
to the strong quark self-energy amplitude A(p) > 1 in the ID model which is known 
to significantly slow the approach to Stefan-Boltzmann thermodynamics 0. Such 
effects should be investigated further with a finite range interaction and the rank-2 
version 



10| of the present separable model might help in that regard. 
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